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Table 2: Best fit parameters for low energy {E < 5 keV) SISO+1 spectral fits (PL — power-law, BB = 
blackbody, ED = absorption edge, GA — gaussian line, BREMS = bremsstrahlung, DISKBB = accretion 

disk multitcmperature blackbody) 



Model 






Eovii (keV) 


TOVII 


Eoviii {keV) 


TOVIII 


E^i {keV) T3 or 


PL + 


BB 




61+"-^^ 

^•"^-0.10 












PL + 


BB - 


1- 1 ED 


71+0 08 


0.752 ±0.016 


0.23 ±0.04 








PL + 


BB - 


1- 2 EDs 


"■^+-0.09 


0.74 ±0.01 


0.36 ±0.04 


0.92 ±0.02 


0.22 ±0.04 




PL + 


BB - 


1- 3 EDs (E)* 


1.1 ±0.6 


n 7QO+0.014 
U. (O0_o.095 


S5+° °i 






0.96 ±0.02 0. 


PL + 


BB - 


1- 2EDs + GA {L)i 




0.734 ±0.013 


S5+° °* 

^■"^0-0.03 


0.89 ±0.02 


31+0 05 


0.94 ±0.02 


PL + 


BREMS + 2ED + GA 




0.737 ±0.013 


0.29 ±0.05 


0.90 ±0.02 


n 99+0.06 
"•^^-0.04 


93+0 07 


PL + 


DISKBB + 2ED + GA 


1 2+1-1 


0.734 ±0.012 


-51 +0 05 


0.90 ±0.02 


0.26 ±0.05 


0.93 ±0.02 


PL + 


PL ^ 


- 2ED + GA 


o 9+0.9 

•J-^-o.s 


+0.014 
^•' + '-0.016 


0.21 ±0.06 


Q 09+0 050 
'-'■^^-0.013 


1 5+°°^ 

^•^"J-0.04 


0.94 ±0.10 



tflux of the soft excess component in the [0.4-2 keV] band expressed in units of 10 erg s ^ cm 
referred to as "Model E" and "Model L" in text respectively 
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Table 3: Best-fit parameters and values for High energy {E > 2.3 keV) SISO+1 and GIS2+3 

simultaneous spectral fits 



Model 


r 


Si {keV) 


ai {eV) 


EWi (eV) 


E2 (keV) 


EW2 {eV) 


E3 {keV) 


EW3 {keV) 




SISO+1 




















PL + IGA 


1.84 ±0.06 


6.3 ±0.2 


5+°-^ 


otcn+170 










55; 


PL + 2GA 


1.85 ±0.07 


6.2 ±0.2 


0.7 ±0.4 


320 ± 190 


(■ 44+0.05 


40 ±30 






54( 


PL + 3GA 


l-^4_oo4 


6.3 ±0.2 


5+° '* 


2Yn+190 


fi 44+0.05 


40 ±30 


r. 4+0.7 


21 ±17 


54: 


GIS2+3 




















PL + IGA 


1 sq+0.03 


6.0 ±0.3 


0.9 ±0.6 












60' 


PL + 2GA 


1.85 ±0.05 


6.3 ±0.2 


Q.^{fixed) 


250 ± 80 






5.5 ±0.7 


< 60 


61' 


PL + 3GA 


1 Qo+0.05 
i.oo_Q 03 


6.41 ± 0.08 


0{fixed) 


160ii° 






5.62 ±0.12 


(50+30 

'JU_20 


60( 
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Table 5: Log of temporal subdivision of the observation, SISl count rate (mean count rate =1.11 cts/s), 
OVII and OVIII absorption edges optical depths and [2-10 keV] photon index in the E model case 



Orbit 


^start ^end (-^s)"'" 


Exposure (s) 


CR 


TOVII 


TOVIII 


r2.3-10 keV 


1 


0-20 


10448 


1.60 


0.28 ±0.08 


0.16 ±0.07 


2.08 ±0.09 


2 


20-60 


15105 


0.90 


0.37 ±0.09 


0.17 ±0.08 


1 68+°°^ 


3 


60-64 


1382 


1.48 


< 0.57 


< 0.35 


2.0 ±0.2 


4 


64-75 


5291 


1.14 


0.48t°;l? 


0.22 ±0.12 


1.76 ±0.09 


5 


75-81.5 


3521 


1.28 


0.23^°:!^ 


0.18 ±0.14 


1 oq-l-0.54 


6 


81.5-88.5 


3884 


0.96 


< 0.36 


< 0.26 


1.82 ±0.11 


7 


88.5-95 


3231 


1.22 


0.30t°;22 


0.16 ±0.16 




8 


95-108 


5281 


0.81 


0.32 ±0.16 


< 0.18 


1.81 ±0.03 


9 


108-115 


2499 


1.31 


0.34 ±0.22 




2.1 ±0.3 


10 


115-140 


10289 


0.66 


63+°- 


< 0.16 


1.6 ±0.4 


11 


140-155 


6646 


1.50 


0.35 ±0.11 


< 0.13 


1 Q1+0-07 

-'-•yj-_o.o6 



'^times are calculated from the start of the observation 
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Table 6: Best-fit parameters for the low energy spectral features in SISO+1 spectra of time-resolved 

plias(^s 



Phase 


Eovii {keV) 


TOVII 


Eoviii {keV) 


TOVIII 


E3 (keV) 


rs or EW3 {eV) 


Xi/ 


Model L 
















High State 


0.726 ±0.012 


0.29 ±0.07 


0.88 ± 0.04 


0.30 ±0.06 




8tl 


1.06 


Low State 


n 7QQ+n.oi3 
u. (oy_o.oo8 


n 44+0-08 


n 90+""'"' 

"•^^-0.02 


n 01 +0.09 
"••J^-0.11 




< 14 


1.06 


10th orbit 


0.722 ±0.017 


0.66 ±0.12 


91+0 02 

"-'•^-'■-0.03 


0.44 ±0.10 


0.94 ± 0.08 


5t 


1.01 


Low State* 


0.743 ± 0.013 


n QQ+0.09 
"•"^^-0.10 


0.87 ±0.04 




90+° °2 
'J.yu_o 05 


16t^o 


0.99 


Model E 
















High State 


0.723 ±0.013 


29+° °'^ 


0.86t 


0.14 ±0.07 


0.96 ±0.03 


0.14 ±0.05 


1.06 


Low State 


n 700+0.015 

U.<00_Q 007 


0.45 ±0.08 


0.86+ 


< 0.13 


0.945tro^ 


^•^-■^-0.05 


1.06 


10th orbit 


r, 799 + 0.017 
U. 1 ZZ_Q 010 


0.66 ±0.12 


0.86t 


< 0.17 


0.961°:°^ 


0-38_o+4 


1.01 


Low State* 


U./4D_Q 013 




0.86t 


< 0.20 


0.94 ± 0.03 


'-'•^■'■-0.08 


0.96 



+ unconstrainted 
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Table 1 : - Good time interval selection criteria 



SISO+1 GIS2+3 

SAA.eq.O 
T_SAA.gt.60 
BR_EARTH.gt.20 
ELV.MIN.gt.5 
T.DY_NT.gt.5 
COR.MIN.gt.6 COR.MIN.gt.7 
SX_PIXLY.lt.400t GZ_Ll.gt.Ot 
tX=0,l, Y=l,3 
tZ=2,3 

SAA = South Atlantic Anomaly (0 = passage excluded) 
T_SAA = Time after passage through SAA (in s) 
BR_EARTH = Angle between pointing direction and Earth's terminator (in degrees) 
ELV_MIN = Angle between pointing direction and Earth's limb (in degrees) 
T_DY_NT = Time after day /night transition (in s) 
COR_MIN = Minimum cut-off rigidity allowed (in GeV/c) 
SX_PIXLY = Maximum number of counts in chip Y of detector SISX 
GX_L1 = GIS2 and GIS3 discriminator 
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Table 7: Best fit parameters for SISO+1 and GIS2+3 simultaneous high energy fits on A02 data quartiles 
with constant direct vs. reflected normalizations ratio (quartile's definition in text, in brackets the [2-10 

keV] luminosity in units 10*^ erg s~^) 



Phase 


r2.3-10 keV 






EWf,:, {eVf 


Ts-lO keV keV 


X7d-0.f. 


F-test 


1st quartile (0.33) 
2nd quartile (0.56) 
3rd quartile (0.70) 
4th quartile (1.13) 


1.57 ±0.08 

1.72 ±0.07 
1.88 ±0.07 


6.28 ±0.14 

6.58 ±0.14 

-0A6 

6.30 ±0.10 


o.o8^^:^,I 

'-'•■'-'^-0.13 
n 4+0-5 


200 ± 90 
270 ± 120 


1.19 ±0.19 

1 n«+ni6 

-■^•^8-0.20 

0.88tn^ 


295/265 
299/294 
301/344 
578/581 


2.4 
0.1 
0.3 
2.5 



^ the quoted F-test is calculated after allowing the J7/27r parameter to float in the reflection model. The 
resulting best-fit values for this parameter are listed in the adjacent column 

and a are referred to the main component of the emission complex; the EW instead is referred to the 
complex as a whole 



6 



Table 8: Best-fit parameters for SISO+1 and GIS2+3 simultaneous high energy fits on A02 data quartiles 
witli refi(X't{xl eompoiK^nt norniahzatioii invers(^ly correlated witli total flux 



Phase 


r2.3-10 keV 






xVd.o.f. 


0/27r (fixed) 


1st quartile 


1.66 ±0.09 


6.0 ±0.5 




288/277 


1.72 


2nd quartile 


1 70+0.22 
' ^-0.10 


6.6 ±0.4 


J. (.+180 
■'■■'^"-110 


290/277 


1.2 


3rd quartile 


1.88 ±0.07 


6.2 ±0.2 


200 ± 90 


322/277 


1 


4th quartile 


2 n^+o-o^ 

Z.UO_Q 06 


5.8 ±0.5 


150 ± 100 


290/277 


0.77 



Hine parameters refers to the main component of the emission complex 
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Table 4: Best-fit parameters for the line emission from a relativistic accretion disk model 



q R^/Rs Ro/Rs 


e 


EW {eV) 


xVd.o.f. 


= 6.4 keV (cold disk) 
-it 28tl 
-3t Wtl 34l~ 
-2 ±3 lot 30t 


25+1 
28tt 
27 ±4 


270 ± 40 
27Qtll 
270 ± 50 


1157/1230 
1157/1230 
1157/1230 


Ea = 6.7 keV (ionized disk) 
-3t 7 ± 3 5000* 
-2.9 ±0.2 8t 5000t 




260 ± 50 


1157/1230 
1157/1230 



t fixed 
t unconstrained 
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Abstract 



We present the results of an extensive analysis of the ASCA A02 observation of the Seyfert 1 galaxy 
NGC4051. The target exhibits broadband [0.5-10 keV] variability by a factor 8 on time scales ^ 10^ s, 
with a typical doubling time ~ 500 s. The spectrum is characterized by a strong emission excess over 
the extrapolated power law at energies E < 1 keV. Absorption edges due to ionized oxygen species 
OVII and OVIII are detected together with an emission-like feature at -E ~ 0.93 keV. The OVII edge 
undergoes significant variability on a timescale as low as 10^ s, whilst no contemporary variability 
of the OVIII feature is detected. Typical variability time scales place constraints on the location and 
the density of the absorbing matter. In the self-consistent hypothesis of a high energy {E > 2.3 keV) 
power law reflected by an infinite plane-parallel cold slab, a photon index change (AF = 0.4) has also 
been observed; a natural explanation can be found in the framework of non- thermal Comptonization 
models. The iron line is redshifted (centroid energy E 6.1 keV) and broad {a > 0.2 keV); multi- 
component structure is suggestive of emission from a relativistic accretion disk; however if the disk is 
not ionized a contribution by a molecular torus or an iron overabundance by a factor ^ 1.5 are required. 



Key words: Galaxies: individual (NGC4051) - Galaxies: Seyfert - Galaxies: X-rays 



1 Introduction 

NGC4051 is a nearby {z — 0.0023) Seyfert 1 galaxy well known as a rapidly variable 
X-ray source. It was extensively observed in the past decade and both its spectral 
features and its variability deserved great attention. An emission excess above the 
extrapolation of a simple power law was discovered by the Einstein (Urry et al. 1989) and 
E'XO.S'^IT' (Lawrence et al. 1985) experiments in the soft X-ray band. Recent analysis of 
two pointed i?(95'74r observations indicated the presence of a strong absorption feature 
at energy ~ 0.85 keV (M^Hardy et al. 1995); it was interpreted as a blending of 
the photoionization absorption edges of OVII and OVlll. Similar features have been 
commonly discovered in the i?05'74r spectra of Seyfert 1 galaxies and explained as due 
to the presence of highly ionized matter surrounding the central source (the so called 
"warm absorber" model); they have been confirmed by several ASCA targets (among 
the others MCG-6-30-15, Fabian et al. 1994, NGC 3783, George, Turner & Nctzcr 1993, 
NGC3227, Ptak et al. 1994, NGC5548 and Mrk841, Otani 1996). M^Hardy et al. (1995) 
calculated self-consistently spectra emitted by matter in photoionization equilibrium 
with the ionizing continuum, assuming a spherical shell geometry and uniform density; 
they yielded a value for the ionization parameter ^ ~ 300 and a column density Ny,arm. ~ 
8 X 10^^ cm~^, with an intrinsic energy index ~ 1.19. 

At higher energy, there was also much evidence of deviation from the simple power 
law trend. Ginga data have been widely analyzed and many different models were 
claimed to be the best explanation for the observed spectra. Matsuoka et al. (1990) sup- 
ported an interpretation in terms of a partially covered power law with solar abundance 
or a reflection model or a thermal bremsstrahlung and a power law, whilst Kunieda et 
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al. (1992) favoured reprocessing by cold blobs. The detection of an iron K-fluorescence 
line by neutral or mildly ionized matter provided further support to the hypothesis of 
nuclear radiation reprocessing. 

NGC4051 was observed by the ASCA sateUite during the Performance Verification 
phase (PV hereafter). The source fell into one of the SISl inter-chip gaps during it. 
The relevant analysis was published by Mihara et al. (1994) (hereafter Paper I) and 
their main conclusions can be summarized as follows. The unprecedented combination 
of energy resolution and broadband coverage allowed to discover two distinct absorption 
edges at energies Ei ~ 0.74 keV and E ~ 0.92 keV; they were interpreted as due to 
highly ionized oxygen species OVII and OVIII, although the energy of the latter was 
blueshifted (~ 6%). That gave the first direct measurement of warm absorbing gas in 
NGC4051. A power law modified by a warm absorber could partly explain the apparent 
soft excess, but the adding of a further thermal component with kT ~ 0.1 keV was 
required to remove it completely. The best-fit underlying power-law photon index was 
~ 1.88. A narrow emission line at energy consistent with K-fluorescence by neutral 
or mildly ionized iron {E ~ 6.45 keV) was also detected with an upper limit on the 
FWHM of - 460 eV and equivalent width EW ~ 170 eV. 

In this paper we report the results of an extensive analysis of a second ASCA obser- 
vation performed during the A02 phase, with a much longer integration time than the 
PV phase. Timing and spectral analysis and short term (< 2 days) variability will be 
mainly addressed and comparison with the PV phase results will be quoted only in such 
a framework. In §2 data reduction methods and criteria are presented. In §3 we will 
discuss briefly the broadband variability and characterize it from a statistical point of 
view. §4 will deal with the spectral analysis of the A02 data, while in §5 the variabihty 
of the main grand spectral features observed will be studied in detail. §6 will be devoted 
to a discussion of the results, whilst some concluding remarks will follow in §7. 

2 Observation and data reduction 

The ASCA sateUite (Tanaka, Inoue & Holt, 1994) observed NGC4051 during the A02 
phase between 6th July (15:28:34 UT) and 9th July 1994 (10:40:27 UT), for a total 
exposure time of ~ 1.55 x 10^ s. The Sohd State Imaging Spectrometers (SIS hereafter) 
were operating in 1 CCD mode and the data were collected in FAINT mode. Dark Frame 
Error (DFE) and echo corrections (Otani & Dotani 1994) were applied and the Charge 
Transfer Inefficiency (CTI) corrected Pulse Invariant SIS data used. Typical DFE values 
at the time when the observation was performed were ~ 1 2 ADU, corresponding to 
only ~ 7 eV, while the energy deficit for transfer was ~ 2 3 x 10~^ per transfer both 
for SISO and SISl chips used, which causes about 1% difference of the energy. Another 
source of gain uncertainties is the Residual Dark frame Distribution (RDD); for 1-CCD 
mode such effect is smaller than 5 eV even in the post-launch calibrations. The uncer- 
tenties on the energy scale that can be associated with such corrections can be therefore 
estimated as < 10 in the lowest energy ASCA SIS band and < 20 eV around 6 keV. 
A systematic error of such order has been taken into account in the forthcoming results. 
The Gas Imaging Spectrometers (CIS hereafter) were operating in PHA mode. Good 
time intervals for the data analysis have been selected according to the criteria in Table 
1. They amount to a total integration time of ~ 6.8 x 10^ s for both detectors. This 
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is more than twice the effective integration time of the PV phase observation. Photons 
were extracted from a circular area ~ 4.7' of radius for the GIS and ~ 3.2' for the 
SIS. Data preparation and selection were performed via the XSELECT package v.1.2; 
spectral fits were carried out using the XSPEC program, v. 8. 5 (Arnaud et al. 1991). 



3 Temporal behaviour 

Broadband light curves for the SISO and GIS2 detectors are displayed in Figure 1 (bin- 
ning time At — 312 s). Strong variability by a factor up to ~ 8 on time scale < 10^ s is 
clear. Two fiarc-likc structures can be recognized between 10^ and 3 x 10^ s from the 
beginning of the observation and in the last 10"^ s. Events of very rapid variability with 
a doubling time of ~ few hundreds seconds arc also relatively common. 

X-ray variability on very short time scales has been reported by many authors. Dou- 
bling time scales of few ~ 10^ s were found by Marshall et al. (1983) with the Imaging 
Proportional Counter on board the Einstein satellite, by Lawrence et al. (1985) with 
the LE and ME EXOSAT detectors and by Matsuoka et al. (1990) in the 1987 Ginga 
observation. In Paper 1, single episodes of flux doubling on time scales of the order of ~ 
few hundreds of seconds were reported. In order to characterize such variability from the 
statistical point of view, an analysis of the lowest-doubling times has been performed on 
the 100 s-binned light curves of all the four detectors. The results for the SISO and G1S2 
are shown in Figure 2. The dNevents/d{A.t) distribution exhibits a flat behaviour with a 
sudden decrease below some typical time scale r. Fit of the distribution function of the 
form A[l — e*^^^*/^-*], where A is a normalization constant equal to the asymptotic mean 
of the distribution, yielded the values: r = 500 ± 300 s for SISO and r 400 ± 300 s for 
GIS2 light curve, with a reduced ~ 1. We therefore conclude that NGC4051 presents 
broadband statistically significant variability on timescale as low as ~ few hundreds of 
seconds. 



4 A02 observation spectral analysis 

Previous X-ray observations of NGC4051 showed a very complicated and variable spec- 
tral structure; in PV phase data a simple power law with photoelectric absorption by 
cold matter with cosmic abundance was too a simple model to explain the various fea- 
tures observed. 

A02 data confirm such outcomes in the light of a better S/N due to the more 
than doubled effective integration time. Background subtraction has been performed 
both with blank sky data files and with spectra obtained from a source-free region in 
the same FOV of the source. The results are fully consistent and the ones obtained 
with the former procedure are shown in the present paper. Spectral analysis has been 
restricted to the bands [0.57-10 keV] for the SIS and [0.7-10 keV] for the GIS detectors 
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in order to avoid low effective area energy cliannels and tlie systematic narrow-band 
spectral feature observed in the SIS spectra around the K neutral oxygen absorption 
edge {E ~ 0.54 keV, Gendreau et al. 1995) that could significantly affect the results 
in the softer band. In Figure 3 the results of a fit with a simple photoelectric absorbed 
power law are shown. The Nh value has been constrained throughout the paper not to 
be lower than the Galactic value of 1.3 x 10^° cm~^ (Dickey & Lockman 1990). The 
reduced is rather poor in both cases {Xsiso+i = 1644/601 d.o.f., XgiS2+3 = 1633/1304 
d.o.f.). The main contributions to the are due to: a) a soft X-ray excess emission 
above the extrapolated power-law below ~ 1 keV; b) a broad emission line feature 
centered at £' ~ 6.1 keV. 

The average source unabsorbed flux - calculated in the hypothesis of the phenomeno- 
logical "L" model (see §3.1) - is 2.6 x 10"^^ erg cm-^ s'^ in the 0.4-2.0 keV band and 
2.4 X 10~^^ erg cm~'^ in the 2.0-10.0 keV band (average on the SIS spectra), cor- 
responding to source-frame luminosities of 0.58 x 10^^ erg and 0.54 x 10^^ erg 
at z = 0.0023. Here and hereafter Hq = 50 km s^^ Mpc~^ and q = 0.5 are sup- 
posed. Such results are consistent with the 2-10 keV flux range observed by previ- 
ous experiments (0.2 5 x 10~^^ erg cm~'^), including the PV ASCA observation 
(2.2 X 10-^^ erg s'^ cm-^). 



4.1 Soft excess 

In Paper I, Mihara et al. (1994) identified absorption edges features due to highly 
ionized oxygen species OVII and OVIII in the SIS NGC4051 spectra of PV phase data. 
The adding of such features to the simple absorbed power-law best-fit model, although 
statistically justified according to the F-test, is not able to remove the soft excess. 
A complete removal of the soft excess via narrow features (absorption edges and/or 
narrow emission lines with Gaussian dispersion a fixed equal to 0) requires 6 lines in 
the energy range 0.57 < < 1.16 keV and no edges with a = 661/594 d.o.f. 
Such explanation, although statistically acceptable, is de facto simply a simulation of 
a continuum soft excess emission due to the limited energy resolution available; thus it 
will not be considered in the following. 

Instead, adding a thermal continuum component {i.e. a blackbody) to the SIS 
spectrum produces a dramatic decrease of the with Ax^ = 876. However, even after 
this improvement narrow absorption features are clearly visible in the spectrum and 
can be modelled either with a couple of absorption edges and an emission line (hereafter 
Model L) or with three absorption edges (hereafter Model E) , with comparable resulting 

(xi = 624.4/598 d.o.f., x| = 624.0/598 d.o.f.). The adding of each of the quoted 
features is justified according to the F-tcst as shown in Table 2, that contains the 
best-fit parameters for the models described. The fits in Table 2 have been performed 
on the spectral channels in the range [0.57-5 keV], in order to avoid contamination 
by the higher energy line structure. Here and hereafter errors are at 90% level of 
confidence for one interesting parameter (Ax^ = 2.71, Lampton. Margon & Boyer 
1976). This underestimates the actual 90% confidence region for multi-parameter fitting 
but nevertheless serves to provide an indication of the curvature of the x^ space near 
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the best-fitting points. Parameters with units of energies are quoted in the source 
reference frame. The energies of the softer two absorption edges {Ei ~ 0.74 keV and 
E2 — 0.87 keV) are fully consistent - within the statistical uncertainties - with the 
K-shell photoionization energies of the oxygen ion species OVII and OVIII and will be 
therefore associated with them in the following. The depth of the OVII edge is almost 
the same in Model L and E (r ~ 0.35) while the OVIII depth is model-dependent, 
ranging from t'^^^ ^ 0.12 to t^^^ ~ 0.31. The narrow spectral features' properties are 
not significantly dependent on the shape of the continuum, as shown in Table 2 where the 
best-fit parameters obtained when the blackbody is substituted by a bremsstrahlung, 
power-law, or multitemperature disk blackbody emission are shown. The temperature 
for optically thick emission is kT = 180 ± 20 eV. The best-fit blackbody parameters 
correspond to a typical emitting region size Rbb ~ 1.3 x 10^° cm, while the emission 
measure for optically thin thermal emission is nlV ^ 2 x 10^^ cm. 

It's worthwhile to note that a standard scattering model {i.e. the superposition 
of power-law spectra with the same indices seen through different absorbing column 
densities) fails to reproduce the observed soft excess {xf — 1-9, with significant residuals 
in the whole E < 2 keV range). 

If the average SIS spectrum absorption features are fit with only two absorption 
edges, as in Paper I, their energies turn to be Ei — 0.738 ± 0.013 keV and E2 — 
0.927 ± 0.017 keV. This outcome is consistent with Paper I; however, according to the 
higher S/N A02 data, E2 is more likely to be a blending of a double structured feature. 
If the model E (L) is applied to the PV phase data, the resulting optical depths are 
Tovii = 0.261°:?? {tovii = 0.431°:°^) and toviii = 0.13 ± 0.10 {toviii = 0.51 ± 0.04). 

Some of these features he within the useable energy bandpass of the GIS detectors. 
The soft excess in the GIS2 average spectrum can be modelled with a single absorption 
edge at E = 0.94 ± 0.05 keV, with r = 0.34 ± 0.15 and P = 1.88l|j:{]^ (F-test = 163); 
but fits with a comparable statistical significance can be achieved by the adding of 
an emission line at E = 0.81 ± 0.03 keV {EW = 70tlo ^V, P = LSGtofe) for GIS2 
and E = O.SOto^ol eV, {EW = 60tlo eV, P = 1.88l[j:[]^) for GIS3 or with a thermal 
continuum. 

4.2 Warm absorber model 

As already stated in §1, absorption features from highly ionized oxygen have been widely 
observed in the RO SAT and ^l^C^l spectra of many Seyfert 1 galaxies and suggest the 
presence of warm matter along the line of sight. Models of absorption by ionized matter 
intervening along the line of sight (the so called "warm absorber") can be calculated 
and compared with the available observational data, provided some simplyifying as- 
sumption on the geometrical and physical properties of the absorbing matter are made. 
We have used the photoionization code Cloudy (Ferland 1991) in order to calculate 
self-consistent total spectra resulting from the passage of ionizing continuum radiation 
through an optically- and geometrically-thin spherical shell. A power-law shape has 
been assumed for the incident spectrum in the whole energy range between 13.6 eV and 
40 keV. A complete characterization of the physical properties of the absorbing matter 
is given once the parameter ^ = L/ueR^ is provided, where L is the luminosity of the 
ionizing continuum, Ug the numerical electron density and R the distance between the 
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central source and the inner edge of the shelL In fact the physical properties of the 
absorber matter depend on the ionizing flux F ca L/ B? = rie and ^ are supposed 
uniform inside the shell, while R — 10^^ cm and L — 10^^ erg s~^. The free parameters 
in such models arc the photon index F, the ionization parameter ^ and the "warm" 
hydrogen column density Ny/. A grid table model has been fitted to the SIS spectrum. 
Figure 4 shows the best-fit data/model ratio, with = 851/604 d.o.f. The major 
contributions to the below ~ 1 keV, where the warm absorber features are 

present that should provide the bulk of information on the ionization state of the matter. 
The best-fit parameters are: F ~ 2.12, log{Nw) ~ 22.27 and ^ ~ 110. A prominent 
emission-like feature is present at ~ 0.90 keV and the OVll edge looks not to be well 
modelled. An acceptable can be obtained either through the adding of a thermal 
component, like a blackbody (x^ = 694/602 d.o.f.) or through the adding of three emis- 
sion lines at energies E = 0.60 ± 0.08 keV , E = 0.69 ± 0.08 keV , E = 0.93 ± 0.07 keV 
(x^ — 660/600 d.o.f.), but in both cases the warm absorber matter best-fit parameters 
are not consistent with the presence of both OVll and OVlll ions [for example, in the 
latter: Ioq^Nh) > 22.50, ^ = 162 ±12]; such solutions, although statistically acceptable, 
have to be regarded only as purely mathematical. Thus the model does not provide a 
satisfactory explanation of the data as the approximations and assumptions appear to 
be too simple. A more detailed discussion on the properties of the absorbing matter 
implied will follow in §5, in the light of the variability properties that it displays. 

4.3 Iron line emission 

Adding a Gaussian emission line to the average spectra of both SIS and GIS detectors 
(in order to model the residuals around the Fe-K region) improves significantly the x^- 
The best-fit parameters for a simple absorbed power-law and Gaussian line model are 
shown in Table 3, where only the channels with £■ > 2.3 keV have been included in the 
fit in order to avoid contamination by the soft excess spectral components. The line 
centroid energy: Esis = 6.3 ± 0.2 keV and Eqis = 6.0 ± 0.3 keV; the equivalent width 
is quite large compared to the typical values observed in Seyfert 1 galaxies {EWsis = 
350II50 eV) when compared to Ginga results. The line has a clear broad structure. 
Power-law photon index and normalization are F = 1.841q;q| and A^^^ = 7.3lo;| x 
10~^ photons keV cm~'^ at 1 keV respectively. 

A detailed analysis of the emission feature reveals a much more complex structure. 
X^ is further improved by adding two narrow {i.e. Gaussian dispersion a held fixed at 
0) Gaussian emission components with energies 

^(2) 6.4 keV and E^^) - 5.5 keV , 
albeit the significance of the third is somewhat marginal in the SIS spectra (Axl/^ = 4, 
Axg75 = 11)- The best- fit parameters for such fits are again summarized in Table 3. 
The presence of the new components affects slightly the centroid energy of the main 
component that decreases by ~ 2%. We will refer to such multi-component structure 
as "emission complex" in the following. 

Iron K-fluorescence fine is considered as a typical signature of the reprocessing of the 
central source radiation by optically thick matter, possibly in the form of an accretion 
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disk around the putative black hole. Such a scenario, investigated by the works of 
Guilbert & Rees (1988) and Lightman & White (1988), predicts a continuum "bump" 
at energies £■ > 10 keV. Since the appearance of the line and the excess continuum 
emission are thought to be strongly correlated, a power-law + reflection by a "cold" 
{i.e. not ionized) accretion disk + emission complex model has been applied both to 
the SIS and GIS data. The inclination of the disk has been fixed to = 30° (sec below), 
but the dependency of the following results on this parameter is completely negligible. 
The only free parameter of the reflected spectrum is then the solid angle subtended 
by the reflecting matter Q/2n. Adding of this new degree of freedom is not justifled 
according to the F-test {Fqis = 1-0, Fsis = 0.1); only a very loose constraint on the 
solid angle (f2/27r < 7) is yielded if the model is applied on the spectra of the four 
detectors simultaneously. In the hypothesis of reflection by a plane-parallel slab (f2/27r 
held flxed to 1), the best-flt parameters for the main emission complex component are: 
= 6.1 ± 0.3 keV, (t(i) = O.Ttg;^ eV, EW^^^ = 220t?^g eV, with an intrinsic photon 
index T = 1.88l'^o;o4. The additional complex components have: i?*^^^ = 6.44 ±0.05 keV, 
EW^"^^ = 60±30,'e(3) = 5.50l[|:i^ keV and EW'^^^ = IQtjl eV respectively. The contour 
plot for the centroid energy vs. the Gaussian dispersion of the main component is shown 
in the upper panel of Figure 5; at 90% level of confldence the Gaussian dispersion is 
> 0.2 keV. 

The total equivalent width of the emission complex is much higher than the average 
observed by Ginga in the same target {{EW) NGCi05i = 140 ± 70 eV, Nandra & Pounds 
1994), where a single "narrow" line was supposed. It is also much higher than the mean 
of the observed EW distribution in the complete Ginga Seyfert 1 sample {{EW)Ginga = 
140 ± 20 eV, Nandra & Pounds 1994). However several Seyfert 1 observed by ASCA 
showed clearly a broad fluorescence iron line: IC4329A (0-99% > 0.15 keV, Mushotzky 
et al. 1995, Cappi et al. 1996), NGC5548 (0^99% > 0.1 keV , Mushotzky et al. 1995), 
MCG-6-30-15 (^90% = 0.16^:^:06 keV , Fabian et al. 1994). Complex and broad emission 
structures are to be expected if the emission originates in an accretion disk. As stated 
by many authors (Fabian et al. 1989, Matt, Perola & Piro 1991), the Doppler shifting, 
gravitational reddening and relativistic beaming are effective in altering the Gaussian 
profile of an emission line and producing a typical doubled-horned structure. Compelling 
evidence of such kind of structure in the ASCA observation of the Seyfert 1 MGC-6-30- 
15 (Tanaka et al. 1995) and of the intermediate Seyfert NGC4151 (Yaqoob et al. 1995) 
have been recently discovered, giving significant support for such a scenario. 

In order to check quantitatively the disk hypothesis, line emission from a relativistic 
accretion disk model by Fabian et al. (1989) has been applied to the [2.3-10 keV] spectra 
of all detectors simultaneously. The free parameters are the inner (i?j) and outer {Rq) 
radii of the line emitting region in units of Schwarzschild radii Rs from the center and 
the inclination angle 9 between the normal to the plane of the disk and the line of sight. 
The energy of the line centroid E^ has been held fixed to 6.4 (cold disk) and 6.7 (ionized 
disk). The emerging profile is also affected by the radial dependence of the emissivity, 
that in such a model is parameterized with a function R^'^. The underlying continuum 
is supposed to be a power-law with reffection by a plane-parallel slab and the disk-line 
inchnation tied to the reflecting slab inchnation. The best-flt parameters are hsted in 
Table 4. The observed line is consistent either with a cold disk at moderate inclination 
{9 ~ 30°) or with an almost face-on ionized disk. Ri is ~ IO-R5, while only in the cold 
disk case Ro can be constrained to a value ~ 30i?5. In the cold case the derived EW is 
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higher than the maximum allowed by presently available models with solar abundance, 
by a factor ~ 1.5. 

We stress that no absorption edge by neutral or ionized iron has been detected in the 
average A02 SIS spectra. Upper hmits on the optical depth for some trial energies have 
been: TpeiiE = l.llkeV) < 0.06, TpeXxiE = 8 keV) < 0.19, TpeXxiviE = 8.80 keV) < 
0.05. 



5 Spectral variability 

NGC4051 is known to be a variable X-ray source. Both EXOSAT (Lawrence et al. 1985) 
and Ginga (Matsuoka et al. 1990) observations showed a hardening of the spectrum for 
decreasing flux; in both cases it was described in terms of decrease of the spectral photon 
index. On the other hand, Kunieda et al. (1992) analyzed successfully the spectrum 
of the May 1988 Ginga observation in terms of reprocessing of the nuclear radiation by 
blobs of matter {Nhiob ~ 10^^-^ cm~^). Significant spectral changes accompanied fiux 
variations by a factor of 2-3 within 1000 s in the softer Ginga band (2.3-6.4 keV), but 
they did not strictly follow the luminosity trend; instead, it was possible to distinguish 
three different phases, characterized by a constant hard {i.e. 8.7-20.9 keV) flux with 
hard spectrum, by a low energy flare with soft spectrum and by a soft spectrum with al- 
most synchronized fiux variability respectively. Such complex spectral changes could be 
successfully explained by changes of the intrinsic luminosity and the blob number along 
the line of sight. Fiore et al. (1992a) re-analyzed all the Ginga NGC4051 observations 
in the framework of two-component models (power-law -|- reprocessing) and affirmed 
a simultaneous variation of both. A 99% level of confidence variability of the iron line 
correlated with that of the underlying continuum down to time scales as short as few 
thousands seconds gave further support to such hypothesis. 

The ASGA A02 observation data provide a lot of information about the variability 
of the source and allow to deliver new insight on this topic. In the following the soft 
excess and intrinsic power-law photon index variability will be mainly addressed, both 
via hardness ratio analysis and via direct comparision of spectral fits performed on 
different intensity phases. 

5.1 Hardness ratio analysis 

The broadband spectrum of the NGC4051 A02 ASGA observation shows fiux-dependent 
changes at both low and high energies. The PHA ratio of the average SISO spectra 
corresponding to a High and Low intensity State is shown in Figure 6. An operative 
definition of the procedure followed to select High and Low State photons will be given 
in §5.2. The High State spectrum is harder for E >2 keV and E < 1 keV. 

A useful tool for characterizing the spectral variabihty of X-ray sources is the study 
of Softness Ratio {SR hereafter). The energy bands have been chosen in such a way 
they are representative of the behaviour of different spectral components; according to 
Netzer, Turner & George (1994) the most suitable choice is the bands: [0.7-1.3 keV] and 
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[2.5-5 keV]. The SR is therefore defined as: 



_ CR [0.7-1.3 kcV] 

CR [2.5-5.0 keV] ^ ^ 

Such choice allows to study the differential variability of the spectral components. Warm 
absorber should predominantly affect the ffux in the softer band. The low energy cut-off 
of the harder band has been chosen in order to avoid contamination by the iron line 
emission complex. 

In Figure 7 the SISO-l-1 light curve in the two bands and the SR light curve are 
shown for a binning time At = 500 s. The for constant hypothesis of the SR light 
curve is Xsr ~ 341.2/61 d.o.f. The SR undergoes a significant variability, getting lower 
for increasing flux and higher for decreasing flux. When the source flux approaches 
its minimum throughout the observation (1.2 x 10^ s after the start), the SR doubles 
within a time scale ~ 10^ s. A decrease of comparable order of magnitude seems to 
be associated with the sudden doubling of flux within ~ few thousands of second after 
~ 10^ s from the start. 

The determination of the correct variability timescale is of the utmost importance in 
order to derive meaningful constraints on the geometry and physical processes involved 
in Active Galactic Nuclei (AGN hereafter). The SR dependency on the total count rate 
has then been studied in detail. In Figure 8 the SR is plotted vs. the total [0.7-5 keV] 
count rate (binning time At = 5000). The fo^^ constant hypothesis is = 259/29 
d.o.f., corresponding to a chance occurrence likelihood P < 0.1%. Also the linear 
correlation coefficient is rather high (r29 = 0.73 ±0.05), again corresponding to a chance 
occurrence likehhood < 0.1%. A careful visual inspection at the SR fight curve in Figure 
7 would suggest the highest variability to be concentrated in the ~ 2 x 10^ s time interval 
when the source fiux is minimum. The same analysis as above has then been repeated 
excluding it from the light curves; however the strong SR vs. CR correlation has been 
confirmed (x^ = 131/24 d.o.f., r24 = 0.64±0.10, P < 0.1% for both parameters). Similar 
results are obtained for binning times At = 500 s (x^ = 320/61 d.o.f., r = 0.73 ± 0.05) 
and At = 25000 s (x^ = 151/6 d.o.f., r = 0.95 ± 0.03, P ^ 0.5%). The binning time 
is obviously limited by the total elapsed time of the observation {T^iaps ~ 1.5 x 10^ s), 
and no significant piece of information can be derived for At > 25000 s. Anyway, there 
is compelling evidence of an increase of the SR when the source flux increases on the 
whole two order of magnitude range between At ~ 10^ s and At ~ 10^ s. That evidence 
can naturally be explained in terms of a lower hard flux or weaker absorption by ionized 
matter when the source is more intense or any combination of the two. In order to 
disentangle the two effects (or others that might be present) a detailed time-resolved 
analysis of the various spectral components is needed. 

5.2 Spectral analysis: soft excess 

The high effective integration time allows a detailed time-resolved analysis of the main 
grand spectral features discovered in the NGC4051 average spectrum. For such purpose, 
the whole observation has been divided into 11 contiguous segments according to the 
following procedure: a At = 100 s light curve has been extracted and all the consecutive 
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bins whose count rate was higher then the total mean count rate have been gathered 
together (contamination not higher than ~ 12% has been allowed). The same procedure 
has been followed for bins whose count rate was lower than the total mean count rate. 
The log of the resulting subdivisions is shown in Table 5. Each of such subdivision will 
be referred to as orbit in the following, although they have no relation with the real 
orbits of the ASCA satellite. The span of each orbit is never less than ~ 7 x lO'^ s 
except in the case of the 3rd orbit whose time span is ~ 4 x 10^ s. The count rate for 
each orbit is in the range 60% — 144% of the total mean count rate. The minimum flux 
episode is included in the 10th orbit, the flare-likc structures in the 2nd and list. 

The best-fit optical depths for OVll and OVlll edges when Model L is applied to 
the combined SISO+1 spectrum of each orbit are shown in 5th and 6th columns of Table 
5. In such fits the energies of the edges are frozen to their best-fit values in the average 
spectral fit due to low S/N. The 10th orbit OVII optical depth is significantly higher than 
the mean. For the sake of a deeper investigation, we have gathered together the orbits 
1+3+5+7+9+11 (hereafter referred to as High State, HS) and the orbits 2+4+6+8+10 
(Low State, LS) and applied Model E and L to the SIS HS and LS spectra. The results 
are summarized in Table 6. The OVII optical depth undergoes a significant (90% level 
of confidence) change by a factor ~ 50%, the edge energies remaining consistent within 
the quoted statistical uncertenties. No significant variation is observed either for energy 
or optical depth of the OVlll edge. A further detailed analysis has been performed 
splitting the 10th orbit from the rest of Low State, in order to check the significance of 
the outcomes summarized in Table 5. We will refer hereafter to the Low State phase 
when the 10th orbit is subtracted as LS*. A 99% level of confidence contour plot for the 
OVII edge energy vs. optical depth when the L model is applied to HS, LS* and 10th 
orbit is plotted in Figure 9. The 10th orbit and HS optical depths are different at more 
than 99% level of confidence, whilst the LS* optical depth is not consistent with the 
other two phases' at only 90% level of confidence. At this latter level, however, it can 
be stated that a smooth trend of decreasing of OVII optical depth with increasing fiux 
is observed on timescales as low as 10^ s, and then its variability is not limited at the 
lowest flux interval. No signiflcant variability within the quoted statistical uncertainties 
is instead detected for the OVIII absorption edge (see Figure 10). 

The spectral shape in the lowest energy ASCA band is very complex; the observed 
variability of the narrow-band features might be due to the interplay between the ionized 
absorption and the continuum soft excess component. It must be remembered however 
that the "L" model comprises self-consistently all the spectral components required to 
remove the soft excess. Moreover, the observed variability patterns are almost inde- 
pendent of the detailed shape of the soft continuum component used. In Figure 11 the 
optical depth vs. energy contour plots are shown when a bremsstrahlung or power-law 
models replace the blackbody. Analogous conclusions can be derived about the edges' 
variability. The continuum models in Figure 11 are likely to be less physically mean- 
ingful explanations of the observed soft excess emission; but the lack of any significant 
dependency of the narrow-band feature results on the underlying continuum spectral 
shape gives strong support to their reality. 
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5.3 Spectral analysis: power— law photon index 

The 7th row of Tabic 5 shows the best-fit value that the power-law photon index 
r assumes when the L model is applied. It is significantly variable (x^ for constant 
hypothesis 28/10 d.o.f.) and appears to be correlated with the total flux (rio = 0.8±0.3). 
If the soft excess is efficiently removed by models "E" or "L" , this should be an indication 
of a real variation of the observed underlying continuum photon index. 

Assessing unambiguously an intrinsic variation of the spectral index is however a 
very difficult task if the ASCA data alone are available. Any variation in the relative 
normalization of the direct and reflected components might produce an apparent change 
in the steepness of the spectrum. If, as the outcomes in §4.3 might suggest, the iron line 
comes from a rclativistic accretion disk, its broad and skewed structure might give a 
significant contribution to the photon emission in the 4-8 keV range. Last, but not least, 
absorption by ionized matter might affect the high energy spectrum also above 2.3 keV. 
The hmited ASCA high energy bandwidth does not give a good statistics above ~ 8 keV 
and it is thus difficult to constrain the spectral shape in the hard tail. Having all such 
problems in mind, we searched for index changes, provided a self-consistent physical 
model with some simplifying assumptions is assumed, in order to reduce the number of 
degrees of freedom; namely: 

i) we assumed a power law spectrum -|- reflection by a plane-parallel inflnite slab -|- 
iron emission line complex. 

ii) The fainter components' parameters have been held flxed to their average best-flt 
value (see §4.3) 

Hi) neither the absorbing medium nor the continuum soft excess component have been 
supposed to affect substantially the spectral shape for £^ > 2 keV. That imphes 
the numerical density of the intervening matter < 10^^ cm~^ and the temperature 
of the blackbody of the order at most few hundreds cV. Both conditions are fulfilled 
in the ASCA data according to the low energy bandpass outcomes (see §6.2) 

Such a model has been fitted to the SIS-I-GIS spectra obtained when the following 
orbits are gathered together: 8+10 (1st quartile), 2+6 (2nd quartile), 4+5+7+9 (3rd 
quartilc), 1+3+11 (4th quartile). The relative normalization between the direct and 
refiected components has been kept constant (r2/27r = 1) or allowed to vary inversely to 
the flux; such two cases simulate an immediate or delayed response of the reprocessing to 
nuclear flux changes respectively. In the latter scenario, the intrinsic width of the main 
emission complex component has been held fixed to its average best-fit value in order 
to further reduce the available degress of freedom. The best-fit parameters are listed in 
Table 7 and 8. T shows the same qualitative trend in both cases. In Figure 12 a plot 
of the r2.3_iofeey vs. the [2-10 keV] luminosity is shown in the constant normahzation 
ratio case. The figure includes also two measures when the same model is applied to the 
spectra of the HS and LS PV phase data (the definition of the intensity phase in PV data 
is analogous to the A02 case). If all the six experimental points are taken into account 
the index is variable (x^ for constant hypothesis 33/5 d.o.f.) and correlated with the 
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[2-10 keV] luminosity {tq = 0.88 ±0.08, the lower limit corresponding to a likelihood of 
casual occurrence P < 6%). Such correlation is valid also if only the four experimental 
points relative to the A02 data are considered (x^ for constant hypothesis 31/3 d.o.f., 
r4 = 0.98 ± 0.04, P < 25%). In Figure 12 it is straightforward to see that the PV 
phase measuraments are not aligned exactly with the same straight line corresponding 
to the A02 measures. A secular pattern of variability could therefore be superimposed 
on the short timescale correlation between spectral index and luminosity. If the PV 
phase experimental points are arbitrarily shifted rightward by an amount ~ 20%, the 
alignment is qualitatively much better and the linear correlation coefficient P vs. L 
turns to be consistent with 1 even for all the six measures (re = 0.97 ± 0.04); we remind 
here that the average [2-10 keV] flux of the A02 observation is just 20% higher then 
the flux of the PV phase data in the same energy range. Hysteresis phenomena in the 
P vs. L relation should be expected if the continuum is produced in a pair dominated 
Comptonized non-thermal plasma (Yaqoob 1992). If we take the extreme points in 
Figure 12 as representative of the low and high state of the source, the pivot point for 
spectral change is Ep ~ 20 keV. 

Next question to be addressed is whether the observed changes in the spectral index 
are really due to changes in the slope of the intrinsic power-law. As stated in section 
4.3, the presence of the iron line is supposed to be correlated to a "bump" of continuum 
emission at energies £■ > 10 keV. Differences in the observed photon index could there- 
fore be due to different amount of reflected flux, with the stronger reflection producing 
a flatter spectrum. In order to check this hypothesis we have flrst calculated the ratio 
between the photon indices in the [2.3-5.0 keV] and [5.0-10.0 keV], the former band 
likely to be unaffected by any reprocessed component; the measure of such ratio for all 
the quartiles is reported in 6th column of Table 7. Although a slight trend of decreasing 
of such ratio with increasing luminosity can be noticed, all the values are consistent with 
1. As a further check, flts have been performed when the Q/2'k parameter in the reflec- 
tion model is left free. The adding of this further degree of freedom is not statistically 
justified for any of the A02 quartiles (see 7th column in Table 7) and in neither the HS 
nor the LS phase in PV data. 

We conclude therefore that the AP ~ 0.4 variation observed in the ASCA data of 
NGC4051 might be interpreted as due to an intrinsic index change in the framework 
of a self-consistent physical model that takes into account all the spectral components 
in the refiecting accretion disk scenario. Such outcomes partly confirms the results of 
Matsuoka et al. (1990). 

6 Discussion 

6.1 Broadband variability 

The outcomes presented contribute to shed new fight on the spectral properties of the 
X-ray emission from AGNs and allow to set constraints on the geometrical and physical 
properties of the processes involved. 

The main source for the huge amount of power emitted by the nuclear region of 
AGN is the conversion of gravitational to radiative energy in matter falling onto a 
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supermassive central accretor. NGC4051 is well known as one of the most variable ex- 
tragalactic sources; ASCA observations revealed flux changes by a factor of 8 within 
typical timescale ~ 10^ s and doubling time ~ few hundreds seconds. The ROSAT 
hght curve in the softer [0.1-2.4 keV] range displayed variations by a factor ~ 10 within 
~ 4 X 10^ s (M'^Hardy al. 1995). Light crossing time arguments yield an upper limit 
for the mass of the central object Ad < 10^ Mq if the size of the emission region is as 
large as 5 Schwarzschild radii. We stress that previous estimates based on the corre- 
lation between the iron line emission strength and the underlying continuum intensity 
- that yielded values by a factor of 5 lower - are not confirmed by the ASCA data, 
since no flux-correlated variability of the iron line could be firmly assessed within the 
2 days A02 observations (cfr. Table 7). The required efficiency in converting accreting 
matter to luminosity (Fabian & Rees 1979) is > 10~^, consistent with accretion onto a 
Schwarzschild black hole. 

6.2 Identification of the third narrow low— energy spectral fea- 
ture 

The A02 spectra have provided further confirmations of the presence of narrow features 
due to absorption by highly ionized matter. Better A02 data S/N ratio has allowed to 
recognize the absorption edge discovered in PV phase data at energy E ~ 0.92 keV as 
probably due to blending of two different features. Two edges at energies E ~ 0.74 keV 
and E ~ 0.86 keV have then been detected in the A02 SIS spectrum, whose energies 
are fully consistent with the K-photoionization energies of the ionized oxygen species 
OVII and OVIII. A third narrow feature is statistically required to model the spectrum, 
but it is not possible to distinguish between a further absorption edge or an emission 
line from a purely statistical point of view. A correct identification of the nature of this 
third feature is of primary importance, also because the OVIII absorption edge optical 
depth is strongly model dependent. We will therefore address this topic in this section. 

Provided this third feature is an absorption edge, it could be a blending of the K- 
photoionization edges of Nelll and/or NelV (937 < E < 971 eV) or a blending of 
L-photoionization edges of FeVII-j-X (930 < E < 1000 eV). If such species arc to be 
located in the same physical region as OVII and OVIII, the ionization structure and 
temperature of the corresponding absorbing matter has to be the same. In the "E" 
model the typical values of oxygen optical dephts are tovii ~ 0.35 and toviii ^ 0.12, 
corresponding to hydrogen equivalent column density of Nh{OVII) ~ 1.6 x 10^^ cm~^ 
and Nh{OVIII) ~ 1.3 x 10^^ cm~'^. The range of ionization parameters that allows an 
abundance ratio [OVII] /[OVII I] ~ 1 for a cloud of gas in photoionization equilibrium 
is relatively narrow (^ = 50 -^ 90 erg cm s~^, Kallmann & McCray 1982), with typical 
plasma temperature ~ 2 x 10^ K. On the other hand, the presence of Nelll/NelV or 
FeVII-=-X species implies ^ not to be higher than 25. Moreover the equivalent hydrogen 
FeVII^X column density required would be ~ 8 x 10^^ cm~^. One could suppose the 
oxygen features to be produced in a physically different region than the third absorption 
feature; such complication looks however unnecessary, although it cannot be ruled out 
a priori. 

A plasma in photoionization equilibrium with an ionizing continuum is expected to 
emit a rich spectrum of X-ray lines as the result of collisional excitation, inner shell 
fluorescence and recombination. Line-like features with typical EW ~ 30 — 100 eV were 
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detected in the Einstein SSS spectra of several Seyfert 1 galaxies (Turner et al. 1991) 
and suggested to originate in a hot (T ~ 10^ K) coUisionally excited plasma. However 
subsequently Netzer (1993) pointed out that photoionized gas might be responsible for 
them. One of the most promising candidates in the range of the observed NGC4051 
line(s) (A 13.0 — 13. 5A) is L-shell iron line-complex which can cover the whole range 
between 0.7 and 1.6 kcV (Band eA, al. 1990, Liedhal et al. 1990). However, iron L- 
shell lines should dominate the spectrum for values of the ionization parameters U = 
Fion/^HC > 10, where Fjo„ is the ionizing continuum flux and c is the speed of light. 
That corresponds to ^ > 100, slightly higher than the one inferred by the oxygen ions 
abundance. Line-like features with energy ~ 900 eV could therefore be likely produced 
by OVn and OVIH continua or NeVH/IX lines. 

6.3 Warm absorber variability 

The presence of absorption edges due to highly ionized oxygen advocates for a distribu- 
tion of matter along the line of sight whose ionization structure is - partly or totally - 
sustained by the radiation emitted by the AGN central source. In such physical condi- 
tions, the state of the matter can be parametrized through the electronic temperature 
T and an ionization parameter like ^ or U. Krohk, McKee & Tarter (1981) showed that 
in conditions of thermal and ionization equilibrium matter should be thermally stable 
only in two phases: a cold (T < 10^ K) phase with lower ionization states and a hot 
(T > 10^ K) phase with higher ionization states. However, a recent paper by Reynolds 
& Fabian (1994), stated the existence of a possible intermediate stable state with typ- 
ical T ~ 10^ K and ^ ~ 10^-^"^, which they proposed to be a good candidate for the 
warm absorber. Such state would be consistent with the co-existence of OVII and OVIH 
species; the range of allowed ionization parameters is relatively narrow [AS/S ~ 0.1, 
where S = (2 x lO'^K/T)^] and that allows to study the warm absorber variability in 
terms of fluctuations around a stable state where both OVI and OIX numerical density 
is negiglible if compared to the prevailing oxygen ions. 

The detected variability of absorption features in the ASCA A02 observation of 
NGC4051 allows to put some constraints on the geometrical and physical properties 
of the absorbing matter. At first we suppose density and ionization parameter of the 
absorbing matter to be uniform (one-zone hypothesis) and the ionization structure 
to be driven predominantly by the photoionization equilibrium. In such hypothesis 
the relevant time scales for each ionized species are the radiative recombination and 
the photoionization time scale to be compared with the typical variability time scale 
tvar ~ 10"^ s. The most natural explanation for the decrease in the OVII edge optical 
depth with increasing flux is the increasing in the amount of photoionization processes 
OVII— >OVIII. The typical ionization timescale can be expressed as {e.g. Otani 1996): 

~ 4 X 10'^ /novii s 

where novii is the numerical hydrogen equivalent density of OVII and (7ph{Eedge) is 
the photoionization cross-section at the absorption edge energy. From the condition 
10^ s > tyar ^ 'tph it is possible to set novii > 4 x 10^ cm~^, that implies a lower limit 
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for the density of the absorbing matter n > 2.1 x 10^'^ cm~^ assuming cosmic oxygen 
abundances, tph can be also expressed as a function of physical parameters of the warm 
absorber: 

t^r' ~ ^ORleL^s s (3) 

where Riq is the distance of the warm absorber from the nuclear source in units of 10^^ cm 
and L43 the X-ray luminosity in units of 10^^ erg s^^. The variability constraint implies 
therefore Riq < 7. If the change of the photoionization opacity produces a change in 
the column density AA'^^^ ~ AR n, it follows AR < 3 x 10^^ cm <^ R. That should be 
suggestive either of geometrically thin or - more likely - of a quite small filling factor 
for the distribution of the absorbing matter. 

A change in the column density of the absorbing OVII should be accompanied by an 
opposite change in the OVIII column density. However, in the ASCA NGC4051 data 
there is no clear evidence of such change. The relevant time scales involved for the 
OVIII balance are the recombination timescale tree for the OVIII— >OVII process and 
the photoionization timescale t^h^^^ ■ The former is expressed by {e.g. Otani 1996): 

tree (4) 

2^0.76 

~ 2 X 10^ 2 X IQ^^/noviii s 

noviii 

In the framework of the "L" model, noviii ~ •^novii and Uee ~ tph^^ ■ The condition 
of ionization equilibrium for OVII is therefore likely to be fulfilled. In the meanwhile, 
the OVIII photoionization timescale is given by: 

t'^r" ~ 1-3 X lO^Vnoy/// s (5) 
according to equation |^. Then t 

ree ~ ^^h^^^ and that could contribute to keep the OVIII 

species constant. 

Similar evidence of such uncorrelated variation of ionized oxygen edges has been 
observed in the ASCA observation of MGC-6-30-15 (Reynolds et al. 1995, Otani et 
al. 1996), where a significant change in the absorption optical depth of OVIII was 
detected with no appreciable change in OVII. A stratification of the absorbing medium 
was claimed as a possible interpretation. A recent paper by Krolik & Kriss (1995) has 
pointed out that there is no reason a priori for imposing either the thermal equilibrium 
or the ionization balance; whether or not they are fulfilled depends strongly on the 
detailed ionization structure, on the mean luminosity and variability features of each 
AGN. The presence of resonance scattering lines can moreover greatly affect both the 
energy and the optical depth of the observed absorption features, leading to an incorrect 
interpretation of absorbing medium properties. Taking that into account can lead, for 
instance, to a dramatically different interpretation of the low-energy spectral variability 
in MGC-6-30-15, suggesting a fixed column density and ionization equlibrium scenario. 
Some caution must then be employed in the interpretation of the NGC4051 quoted 
results and only high resolution spectroscopy in the soft X-ray band will be able to solve 
all the puzzling questions still open. 



16 



6.4 Photon index variability 

ASCA detected a significant (AF ~ 0.4) and robust change of tfie pfioton index for 
a factor of 4 flux variation. A still higher variability (AF ~ 0.7) had been found in 
Ginga data (Matsuoka et al. 1989), although a possible explanation in terms of variable 
absorption could account for the data equally well (Fiore et al. 1992b). The broadband 
ASCA coverage and the improved high energy resolution have allowed to conclude that 
at least part of such variability can be explained in terms of intrinsic photon index 
variation. 

An analogous change in the power-law slope was found by Leighly et al. (1996) 
in the ASCA observation of the Narrow Line Seyfert Galaxy Mrk766, that showed a 
AF = 0.4 for a flux variation by a factor of ~ 2 on typical timescales ~ few thousands 
of seconds. A renewed interest has thus recently arisen about such topic. 

The most successful model for the continuum production in AGN is the Comptoniza- 
tion of UV or soft X-ray seed photons, which probably originate in an accretion disk, by 
a distribution of relativistic electrons. Depending on the nature of the injected electron 
distribution such models can be thermal or nonthermal. A spherical geometry for the 
interaction region is assumed with radius R and luminosity L, but a fairly independent 
description by the details of the involved geometry can be achieved if the dimensionless 
luminosity to size ratio, or compactness parameter I is introduced (Cavalicrc & Morrison 
1980, Guilbert, Fabian & Rees 1983), defined as: £ = where (Jt is the Thompson 

scattering cross section. A further ingredient in such models is the 7-ray absorption via 
electron-positron pair production; the optical depth of such a process is > 1 for £ > 10 
(Svensson 1994). Provided the constraint on R from the lowest doubling time scale 
R < 1.5 X 10^^ cm holds and the observed 2-10 keV luminosity ~ 10'^^ erg s^^, i'x > 2 
and ^f. of the injected electron distribution can be higher by an order of magnitude if 
NGC4051 reproduces the multifrequency spectral behavior of NGC4151 (Yaqoob 1990). 
The condition for pair reprocessing could then also be achieved in NGC4051. 

The relation between power-law spectral photon index and physical parameters in 
the simple case of an unsaturated thermal Comptonized model has been derived analyt- 
ically by Rybicki & Lightman (1979). F depends inversely on the Compton parameter 
y = T^g. We could therefore suppose the change in the observed [2-10 keV] spectral 
index to be due to the gradual increase of kT with decreasing flux. However, if the 
soft photon input is constant, the pivot point should yield an estimate of the energy 
of the soft photons distribution, which is not compatible with the observed value of 
Ep ~ 20 keV (see §5.3). Moreover, recent calculations by Chakrabarti & Titarchuk 
(1995) have shown that spectral photon index as low as 1.5 can be only achieved for a 
small part of the parameter space with a thermal distribution of nonrelativistic electrons 
and AF = 0.5 requires variations of the accretion rate of two order or magnitude. Such 
kind of dramatic flux changes has never been detected in NGC4051, although they are 
not rare in Narrow Line Seyfert Galaxies. 

A recent paper by Ghisellini & Haardt (1994) demonstrated that it is possible to 
set a one-to-one correspondence between F, T, ie and ig, provided that the thermal 
plasma is pair-dominated. The starvation ratio ie/^s is fixed once the photon index 
is known and the ig spans more than 3 order of magnitude with the variation of the 
temperature of plasma in pair equilibrium. Numerical simulations shows that AF ~ 0.4 
requires a variation of the starvation ratio ~ 40; that implies at least a factor of 10 of 
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variation for the soft-energy seed photon intensity. These numbers are in fairly good 
agreement with the observed ROSAT variabihty (M'^Hardy et. al. 1995); the upper 
hmit on the central black hole mass is two order of magnitude higher than that derived 
from the broadband ASCA variability (see § 6.1) if the ROSAT timB scale corresponds 
to the dynamical time scale of the innermost allowed Schwarzschild orbit. However the 
maximum dynamical range of the soft excess continuum component allowed by ASCA 
data is only ~ 2 (1.0 x 10"^^ < < 2.0 x 10-11 erg s'^ cm-^ if the "L" model is 
applied to the spectra of the quartiles phases). Therefore no evidence of the required 
dramatic flux variability of the soft input component arises from the ASCA data; that 
casts some doubt on the validity of such scenario for NGC4051. 

An alternative explanation is provided by the non thermal Comptonization mod- 
els. With the most simplifying assumptions, soft photons with a blackbody distribution 
peaked at energy Xs = 2.8/cT/meC^ interact with a monoenergetic distribution of rel- 
ativistic electrons with Lorentz factor 70. The compactness parameters ie and ig are 
defined as above. A steady distribution of primary particles produces a spectral photon 
index F 1.5; if every scattered photon makes a pair with equidistributed energies, 
the emerging spectrum of the second generation will have F ~ 1.75. If pair saturation 
is achieved through successive pair cascades, the value of the slope tends to F ~ 2, 
which is consistent with the observed results and could explain the "saturation" effect 
shown by the F2-10 kev vs. L relation in the present analysis. Yaqoob (1992) studied 
in detail the physical condition for the photon index variability observed in the Seyfert 
1.5 galaxy NGC4151 to be achieved. The observed index-fiux correlation in NGC4051 
can be reproduced by a variable Lorentz factor and a constant electron injection rate, 
with correlated input UV variations in order to leave is/^e constant. If the Lorentz 



Similar conclusions were derived in a recent paper by Torricelli-Ciamponi & Cour- 
voisier (1995), but in a slightly different scenario. In their numerical model, both the 
thermal UV emission and the inverse Compton reprocessing take place in the same hot 
plasma cloud. It need not to be small in size nor optically thick to pair production and 
the UV reprocessed photon luminosity is only a fraction of the observed. If the observed 
X-ray continuum soft component is the hard tail of an observed UV, there is therefore 
no need for any correlated variability with the high energy (i.e. 2-10 keV) flux. The 
authors applied their numerical model to the F2-10 keV vs. F relation in Matsuoka et. 
al. (1990) and derived a characteristic size of the cloud ~ 10^^ cm and a UV spectrum 
kT ~ 30 eV. Comparable results might be expected if the fit were to include also the 
outlined ASCA results. But the idea that the thermal and non-thermal radiation to be 
emitted in the same physical region is not strongly supported by the current analysis, 
regardless of the optical properties of the emitting plasma. If we suppose an optically 
thin cloud with a numerical density from the Torricelli-Ciamponi & Courvoisier model 
(n ~ 7 X 10^° cm~^), its typical size is ~ 5 x 10^^ cm; that implies Mbh ~ lO^M©, which 
is one order of magnitude higher than the upper limit imposed by the X-ray variability 
pattern (see §6.1). On the other hand, a shght trend of correlated increase of the thermal 
emission with flux on timescales as low as several 10^ s might suggest to locate an opti- 
cally thick emitting region at ~ 10'^ Rs- In such a case, the characteristic temperature 
T ~ 180 eV is somewhat higher than expected from the innermost region of an optically 



factor varies by a factor of ~ 2 it is constrained to lie between 
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thick disk, but it is fully consistent with the temperature required in a warm absorbed 
power law + blackbody fit of a combined ROSAT and Ginga spectrum (Pounds et al. 
1994). Such outcome might be connected with the low luminosity and black hole mass, 
although it cannot rule out the effect of Comptonization on the inner disk emission. 
Moreover it must be taken into account that the ASCA low energy limited bandwidth 
alone is not the most suitable for the exact determination of the temperature of any 
UV-soft X-ray component and the measured temperature is likely to be an upper limit 
of the physical one. 

We recall that an origin of the soft excess from a frequency-indipendent scattering 
is ruled out by present analysis (see § 4.1). 

The interpretative scenario is however complex and no unambiguous interpretation 
can be claimed as a unique interpretation for the observational data. The results from 
future broadband missions (hke XTE or SAX) will allow to determine with improved 
accuracy the parameters of all the spectral components involved and to constrain the 
possible models and their degrees of freedom. 

6.5 Iron line 

The ASCA A02 observation of NGC4051 detected a clearly broad (cr > 0.2 keV at 
90% level of confidence) and complex emission line, that can be interpreted as a K- 
fiuorescence iron line whose profile has been modified by gravitational and Doppler 
shifts. The profile and centroid of the line allow to derive quantitative constraints on 
the physical and geometrical properties of the emitting region. The emission complex 
observed can be characterized, according to the equations (1) to (3) in Matt et al. 
(1992a): E^a = 6.09 ± 0.03 keV, = 418 ± 16 eV and EWa = 320 ± 30 eV (the 
results arc the weighted mean of the parameters obtained applying the above cquatione 
to each detector continuum-subtracted spectrum). The emission is therefore: a) sig- 
nificantly redshifted (~ 5%) in comparision to the expected energy for K-fluorescence 
emission line from neutral or mildly ionized iron; b) broad; c) particularly intense, be- 
ing the average equivalent width of similar features observed in other Scyfcrt 1 galaxies 
~ 150 eV (Nandra & Pounds 1994). Tanaka et al. (1995) reported the discovery of an 
analogous complex emission feature in the 200,000 ks spectrum of the Seyfert 1 galaxy 
MGC-6-30-15, giving the first conclusive evidence of the effects due to the presence of 
a supermassive black hole in the proximity of the line radiation field. Explanations 
alternative to the relativistic accretion disk seem to be excluded for that object (Fabian 
et al. 1995). Nonetheless the correct interpretation is far from clear and some problems 
are left unsolved even if we consider only the NGC4051 results. Matt et al. (1992) 
showed the EW of a line emerging from an accretion disk can hardly exceed ~ 180 eV, 
a value that is inconsistent with the observed EW ~ 270 eV. An overabundance of iron 
of such order of magnitude could be a straightforward explanation for such discrepancy. 
Moreover, higher equivalent widths can be achieved if the matter in the accretion disk 
is highly ionized (Zycki & Czerny 1994) due to the increased fluorescence yield and the 
reduced photoelectric opacity; EW ~ 300 eV can be achieved if the mean ionization 
state of iron is ~ FeXXIII. In such case the rest energy of the emitted line is expected 
to lie around ~ 6.7 keV. 

The quoted results allow to put some signiflcant constraints on the iron line emission 
region. Following the fully relativistic treatment of Matt et al. (1992a, 1992b), the 
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observed centroid energy and intrinsic width of the hne profile require an inner radius 
~ 10 Rs, with a comparable size of the nuclear source. The outer radius Rg cannot 
be tightly constrained, and might be as high as ~ 10^ Rs if Ri could shrink to ~ QRs- 
In case of a ionized disk, EW ~ 300 eV can be obtained in the hypothesis of an 
extended source around a ~ lO^M© black hole, provided a typical size of the source 
~ 50 Rs and a sub-Eddington luminosity {L/LEdd — 0.5) are achieved. In both cases 
the location of the innermost emitting region would lie in the range 1 -j-5 x 10^^ cm. Any 
response of the line to the variation of the nuclear continuum should lag at most ~ few 
hundreds of seconds. Incidentally, we remark that we have not detected any evidence of 
line variability in the NGC4051 observations, while recent measures claiming correlated 
iron-linc/continuum variability in Seyfert 1 galaxies have recently appeared in literature 
(Iwasawa et al. 1996, Yaqoob et al. 1996). 

Alternatively, a contribution to the line emission from a molecular torus surrounding 
the nucleus has been predicted by some authors (Krolik & Kallmann 1987, Ghisellini, 
Haardt & Matt 1994, Krolik, Madau & Zycki, 1994) and claimed as a possible expla- 
nation for the narrow line observed in the Seyfert 1 galaxy NGC7469 (Guainazzi et al. 
1994, Leighly et al. 1995a). The line should be in such case narrow and centred at 
the nominal K-fluorescence energy for neutral iron [i.e. ~ 6.40 keV). This explanation 
would fit well the multi-component structure of the iron line, that requires a narrow 
Gaussian line at the nominal Ko-fluorescence energy superimposed to the relativistic 
double-horned profile, with EW ~ 60 eV . Such an intensity implies torus column den- 
sity Nh ~ few 10^'^ cmr'^, although much higher densities cannot be ruled out, since the 
EW value is only loosely constrained. 

The discussion above is obviously affected by the choice of the self-consistent contin- 
uum and line model outlined. A confirmation of the sketched scenario with the upcoming 
broadband X-ray missions [i.e. XTE and SAX) is therefore strongly needed in order 
to achieve a deeper understanding of the physical process involved and of the relevant 
geometry. 

7 Conclusions 

The higher the resolution capabilities of the X-ray experiments, the more interesting 

and complex the emission of AGNs appears. The Seyfert 1 galaxy NGC4051 proves 
itself to be an outstanding target, that shares most of the peculiar characteristics of the 
objects of its class, and namely: 

a) strong variability on time scales as short as few hundreds seconds. 

b) Evidence of multicomponent X-ray spectrum, with an apparent soft excess emis- 
sion superimposed on the higher energy power-law. 

c) Detection of absorption edges from ionized oxygen species and emission lines that 
represent clear evidence of warm absorbing matter along the line of sight. 

d) Variability constraints allow to locate the warm absorbing matter at a distance 
from the nuclear source < 10^'' cm with a numerical density n > 2 x 10^° cm~^. 

e) the power-law spectral photon index is intrinsically variable (if reflection by a 
plane-parallel infinite disk is self-consistently supposed, AF ~ 0.4) 
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f) The iron line emission is redshifted, broad and very intense. The outcomes can be 
explained in terms of emission from a relativistic accretion disk, either cold with 
moderate inchnation {9 ~ 25°) or ionized and almost face-on. In the former case, 
the observed EW is higher by a factor ~ 1.5 2 than the maximum allowed by 
presently available solar-abundance relativistic models. The puzzle of EW excess 
might be solved: i): postulating an ad hoc iron overabundance by a factor ~ 1.5; 
ii): if the reflection is from an ionized face-on disk; iii): if a contribution to the 
iron emission complex comes from the putative molecular torus surrounding the 
nucleus. 

The evaluation of the above outlined results must take into account the peculiar 
properties of NGC4051, which is a low-mass {M^h < lO'' Mq), low-luminosity {Lx ~ 
10^^ erg s~^) and strongly X-ray variable (e-folding time ~ few 10^ s) Seyfert 1. The 
mounting evidence for a positive correlation between the 2-10 keV flux and photon index 
point to an interpretation of the nuclear emission in terms of inverse Compton repro- 
cessing of UV seed photons by a non-thermal distribution of relativistic electrons; the 
injected particles are likely to be generated in the immediate proximity of the putative 
black hole, while the exact location of the seed UV emitting region is still controversial. 

NGC4051 exhibits all the main features often seen in Seyfert 1 X-ray spectra. The 
many peculiar features of this source (variable warm absorber, complex iron line struc- 
ture, variable nonthermal continuum slope) allow a deeper insight on the general prop- 
erties of physical processes in AGNs. 

Although a further effort is needed for a full understanding of the complex inter- 
play among the different spectral components and their dynamics (the contribution 
of broadband experiments will be of the primary importance to achieve such a goal), 
some interpretative milestones can be set: a non-thermal continuum, whose production 
mechanism requires an efficient and continuous injection of a non-thermal distribution 
of relativistic electrons, Compton reprocessing of the nuclear radiation, a complex ab- 
sorber with an intrinsic warm component within the inner border of the BLR are the 
common elements any unified scenario must cope with. 
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Figure captions 



Figure 1. Broadband SISO and GIS2 light curves with binning time At — 312.5 s. 
Time is expressed in seconds starting from the beginning of the observation. 

Figure 2. Distribution function of the lowest doubling times for the broad band curves 
of Figure 1 (At = 100 s). The distribution functions have been fit with an empirical 
equation of the form: A[l — ea;p(— At/r)]. Best-fit values are: tsiso = 500 ± 300 s and 
TG752 = 400 ± 300 s. 

Figure 3. Broadband spectra {upper panels) and residuals in units of a {lower panels) 
when a simple power-law with cold absorption is applied. The energy channels have 
been rebinned in order to have a S/N ratio > 25. Soft excess below ~ 1 keV and a 
broad emission line feature around ~ 6.1 keV are clearly visible. 

Figure 4. Data/model ratio when a CLOUDY "warm absorber" model is applied to the 
SISO+l data. Bad modelling of the OVII edge and a prominent emission-like feature 
at energy E ~ 0.95 keV are clearly seen. 

Figure 5. Data/model ratio for a simple power-law model apphed to the SISO and 
SISl spectrum in the energy range [3-10 keV]. The energy channel have been rebinned 
in order to have a S/N ratio > 15. In the upper panel a centroid energy vs. Gaussian 
dispersion contour plot is shown for the emission complex main component when the 
power-law -|- fiat refiection -|- emission complex model is applied to the spectra of all 
detectors simultaneously. The outermost contour corresponds to 90% level of confidence 
for two interesting parameters (Ax^ = 4.61). At such level of confidence the line is broad 
(cr > 0.2 keV). 

Figure 6. PHA ratio between the HS and LS phase SISO spectra (definition of intensity 
phases is given in text). 

Figure 7. [2.5-5.0 keV] and [0.7-1.3 keV] bands and SR SISO+l light curves for a bin- 
ning time At = 500 s. 

Figure 8. SR vs. [0.7-5 keV] count rate for the SISO-l-1 light curve with binning time 
At = 5000 s. The reduced for constant hypothesis is 8.9 and the linear correlation 
coefficient for 29 degrees of freedom r29 = 0.71 ± 0.07. 

Figure 9. OVII edge energy vs. optical depth contour plot for High State (HS), Low 
State subtracted of the 10th orbit (LS*) and 10th orbit. The 10th orbit coincides with 
the ~ 2 X 10'^ s interval when the source fiux is minimum. The outermost contour 
corresponds to the 99% level of confidence for two interesting parameters (Ax^ = 9.21). 
The definition of the intensity phases is given in text. 

Figure 10. The same as Figure 10 for OVIII. 

Figure 11. The same as Figure 9 when a model similar to Model "L" is applied but a 
bremsstrahlung {left panel) or power-law {right panel) emission model is substituted for 
the blackbody. 

Figure 12. Photon index vs. [2-10 keV] luminosity for both A02 and PV phase data 
(indicated by an arrow). 
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